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Mg 2 + -Fe 2+ ORDER- DISORDER AND THE THERMODYNAMICS 
OF THE ORTHOPYROXENE, (Mg,Fe)Si0 3 CRYSTALLINE SOLUTION 

Surendra K. Saxena 
Subrata Ghose 
Planetology Branch 

ABSTRACT 

Activity-composition relation obtained from the intra-crystalline Mg 2 + -Fe 2 + 
distribution data for orthopyroxene, (Mg,Fe)Si0 3 at 500, 600, 700 and 800°C are 
presented. The activity coefficients are expressed as a polynomial in the mole 
fraction X Fe according to Redlich and Kister Equations. The excess free energy 
of mixing is similarly expressed using Guggenheim's Equation: 

" ^Fe ^Mg {^0 + ^1 (^Fe “ ^Mg) + ^2 (^Fe ” ^Mg ) } 

The two constants A x and A 2 decrease systematically with increasing tempera- 
ture. A 0 increases from 800 to 600°C but decreases again at 500°C. The ortho- 
pyroxene solution is asymmetric and the asymmetry increases with decreasing 
temperature. Theoretically orthopyroxene on the Fe-rich side would begin to 
unmix into two coexisting pryoxenes in the vicinity of 500° C. 
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Mg 2+ -Fe 2 + ORDER- DISORDER AND THE THERMODYNAMICS 
OF THE ORTHOPYROXENE, (Mg,Fe)S10 3 CRYSTALLINE SOLUTION 

INTRODUCTION 

Recently Saxena and Ghose (1970) presented the Mossbauer resonance data 
on the distribution of Fe 2+ and Mg 2+ between Ml and M2 sites in orthopyroxenes 
heated at 500, 600, 700 and 800°C. The distribution data fit the "simple mixture" 
model of Guggenheim (1967) for the individual sites. In this paper the informa- 
tion on the solution behavior of I e 2+ and Mg 2+ on each site will be used to arrive 
at the thermodynamic properties of mixing of Fe 2+ and Mg 2+ in the orthopyroxene 
crystalline solution as a whole. The changes in the silicate framework as a 
function of temperature and composition are not included in these calculations. 
Therefore these thermodynamic parameters are not necessarily the same as 
those obtained by calorimetric experiments on the solution of ferrosilite and 
enstatite, but they would be comparable in magnitude. 

ACTIVITY- COMPOSITION RELATION 

The partial activity coefficient of the Fe component with reference to a site in 
the orthopyroxene (Mg,Fe) 2 Si 2 O fi crystalline phase is defined as 

v M1 

f r V = i ir (i - V.‘) a> 
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where /' F M e 1 is the partial activity coefficient, an adjustable energy constant 
and the atomic fraction Fe 2+ /(Fe 2+ + Mg 2+ ) on the Ml site. The partial 
activity of Fe 2+ on the Ml site is 


a 


Ml 

Fe 



(2) 


Similarly for the M2 site. 


a 


M2 
F , 



(3) 


We need a simple model of the macroscopic system, which permits the 
derivation of the activity of Fe component in a crystalline orthopyroxene, (Mg, 
Fe)Si0 3 phase from the partial activities of the component on individual sites 
Ml and M2. On the basis of statistical thermodynamics, such a model is: 


a 


OPX 

Fe 


7 Ml 
Fe 



(4) 


Mueller (1962) used such a relation while discussing a solution model for cum- 
mingtonite. Banno and Matsui (1967) and Thompson (1970) have also used such 
a relationship. 

Figures 1 and 2 show the plots of partial activity of Fe 2+ against mole 
fraction Fe 2 + /(Fe 2+ + Mg 2+ ) on the individual sites. Figures 3, a, b, c and d show 
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the activity-composition relations in orthopyroxene (Fe,Mg)SiO a at various tem- 
peratures. It is assumed that the silicate framework does not change significantly 
and its contribution to the change in activity with temperature may be neglected. 
Therefore a F ° e px is similar to Ojr^/io in these results. Table 1 shows the 
activity-composition data at different temperatures. The relations presented in 
these figures obey the Gibbs-Duhem equation: 





+ x, 


Mg 


d In f M \ 

**Fe I 


PT 
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( 5 ) 


ANALYTICAL EXPRESSIONS FOR THE ACTIVITY-COEFFICIENTS 

It is convenient to express the activity coefficients as a polynomial in If F ° px . 
The following Redlich and Kister equations (see King, 1969) are used for this 
purpose: 


log f, 


Fe 


V. { B - c (3 jt f . - *„)} * D(v - jr„) (5 X " . x M<) m 


f *. = V, {b - C - *■«.)} * ■>(*,. -*■»,) (SJf,. -Jt„) „) 


3 



A least squares curve Is fitted to the data on activity coefficients at each 


temperature using the equation: 

loe (w = b ( j f. - *..) * c ( 6 **' x «. - >) 


+ H*r. - *..) ('- 8 *»,) W 

The values of B f C and D are listed in Table 2. 

Figure 4 shows the plots of log f against I. The curves are based on the 
expressions (6) and (7) and represent the experimental data well, except for 
tnose of the 500°C. 

EXCESS THERMODYNAMIC FUNCTIONS 

Guggenheim (1937) expressed the excess free energy of mixing G E for a 
binary system as: 

CE = *b{ A o * A .(*. - M * M*» - *,)’ ♦•"} (9) 

where A 0 , A, and Aj are related to B, C, and D as A 0 = 2.303 RTE, A 1 = 2.303 RTC, 
and A 2 = 2.303 RTD respectively. The values of A 0 , A t and A 2 are listed in 
Table 2. 
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Figure 5 shows a plot of G E against * Fe at various temperatures. Note that 


G E may also be calculated from the experimental data by using the relation: 


G‘ -- X„ ST to (" ♦ X„ t BT to f,, (10 ) 

G E calculated for the different samples using (10) differs somewhat depending 
upon the fit of the analytical expressions (6) and (7) to the data on individual 
isotherms as shown in Figure 7. 

Excess entropy of mixing, S E and excess enthalpy of mixing, H E can be 
obtained by using the general thermodynamic formulae: 


S E 




d inf 


Fe 


d in f. 


it* * Tf“) -*(*,. * f„ * r„) 


(id 


»* * - RT 'U,.— rr~ * '..-rr 1 ) < 12 > 
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GENERAL CHARACTER OF THE ORTHOPYROXENE, (Mg,Fe)Si0 3 
CRYSTALLINE SOLUTION 

The Asymmetric Solution 

The activity-composition diagrams (Figures 3, a to d) and the plot of G z 
against A' Pe in Figure 5 show that the orthopyroxene solution is asymmetric. 

The asymmetry increases with decreasing temperature. This relation can be 
examined quantitatively by studying the parameters A 0 , Aj and A 3 as a function 
of temperature. Tli8 expression for is: 

^ Fe + ^1 (^Fe ~ ^Mg) + ^2 (*Fe “ ^llg) J (13) 

When the odd terms in the above expression vanish, the solution becomes sym- 
metric. If A 2 and other higher terms are also zero, we obtain the model of 
simple mixture or the "regular" solution with A 0 as the "interchange energy" V, 
The asymmetry in the crystalline solution is then directly a function of the 
parameter A,. 

Figure 6 shows a plot of A 0 , A } and A 2 against temperature. A 0 decreases 
systematically with temperature in the range of 600 to 800*C. A, and A 2 vary 
regularly in the entire range of experimental temperatures. The trend of varia- 
tion for Aj indicates that the asymmetry in the solution increases very rapidly 
with decreasing temperature. On the other hand, with increasing temperature 
A l approaches a constant value of approximately 200 cal/mole. This means 
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that even when A 0 and Aj become zero (approximately at 900° C), that is, when 
the solution is nearly ideal, some asymmetry would still persist. Note that in the 
present case, a positive value of A x means asymmetry of the solution on the Fe 
rich side. 

Free Energy of Mixing and Critical Temperature 

The activity-composition diagram at 500*C shows a tendency for the curve to 
flatten as X ft approaches 0.70. This is related to the asymmetry of the solution. 
It appears that if A, Increases beyond a critical value, unmixing may occur in the 
solution. The critical temperature should lie between 400 and 500° C. Combin- 
ing the analytical expression (9) for G E with the formulae for free energy of 
mixing for ideal solutions, the total ree energy of mixing G M may be calculated 
from the relation: 


G" = G E + RT X rt 4nX rt + RTX U% RTtnX^ (14) 

The values of A 0 , Aj and A 2 at 400 # C obtained by extrapolation are 680, 
1500 and 980 cal/mole respectively. The values of G u at various tempera tureo 
are plotted in Figure 7. The curves tend to straighten out in the Fe rich region, 
especially at lower temperatures. The curve at 40C J C shows a compositional 
range (X ft = 0.75 to 0.90) where the orthopyroxene «v, old unmix into two co- 
existing solutions one somewhat richer in Fe than the other. The critical tem- 
perature when unmixing begins may be close to 500°C. 
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Table 1 


Activity-composition Data on Orthopyroxene, (Mg,Fe)Si0 3 


Y 

500 # C 

600°C 

700®C 

800°C 

Fe 


°Fe 



0.181 

0.236 

0.262 

0.245 

0.216 

0.381 

0.499 

0.492 

0.463 

0.432 

0.455 

0.617 

0.581 

0.549 

0.509 

0.500 

0.694 

0.625 


0.556 

0.580 

0.798 

0.704 

0.668 

- 

0.620 

0.825 

0.738 

0.702 

0.668 

0.720 

0.850 

0.802 

- 

0.755 

0.760 

- 

0.823 

- 

- 

0.860 

- 

0.882 

0.877 

0.872 

v OPX 


°*U 



0.819 

0.838 

0.842 

0.837 

0.828 

0.619 

0.638 

0.672 

0.661 

0.639 

0.545 

0.544 

0.598 

0.593 

0.569 

0.500 

0.472 

0.545 

- 

0.526 

0.420 

0.391 

0.459 

0.473 

- 

0.380 

0.378 

0.439 

0.435 

0.404 

0.270 

0.362 

0.366 

- 

0.321 

0.240 

- 

0.340 

- 

- 

0.140 

- 

0.267 

0.223 

0.186 
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Table 2 


Values of the Constants B, C, D, A 0 , A, and A a for the Analytical 


Expressions for Log f and G E . 


T°C 

B 

C 

D 

A 0 

cal/mole 

Ai 

a 2 

500 

.260 

.300 

.250 

919 

1061 

884 

600 

.280 

.118 

.190 

1118 

469 

759 

700 

.212 

.060 

.130 

944 

267 

579 

800 

.133 

.048 

.067 

651 

236 

329 


11 





Figure 
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Figure 3a. Activity-composition relation in orthopyroxene, (Mg,F«)SiO 

o ^ 

at 500 C. Smooth lines are drawn through the points. 




composition relation in orthopyroxene, (Mg,Fe)SiO. 
Smooth lines are drawn through the points. 








Activity- compos it ion relation in orthopyroxene, (Mg,Fe)SiO. 
o 

at 700 C. Smooth lines are drawn through the points. 
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Figure 4a. 


Plot of log f again* t at 500°C. Open circle* - log f^. 

Filled circle* - log f^. The leaet squares curves are based 
on the analytical expression with three constants B, C and D 
due to Redlich and Kister (King, 1969). 
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Figure 4b. 
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Ot of log f against ^ at 600°C. 
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Figure 5. Plot of the excess free energy G obtained by using Guggenheim's 
(1937) expression with three constants Aq, and A These 
constants are calculated from the three constants B, C and D and 
are listed in Table 2. 
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Figure 7. Plot of the free energy of mixing, G** against Xp e « The solid 
curves are drawn using Guggenheim's analytical expression for 
G E . The open circles represent G** for the samples calculated 
directly by using the equation: G M = Xy^RTln^ + X^RTln^ 

The dashed curve is based on extrapolated values of A^, , 

and A 2 to 400°C and shows that orthopyroxene with Xy g between 
0.75 to 0.90 would unmix into two co-existing pyroxenes. 



